Renormalization group study of a two-valley system with spin-splitting 
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Renormalization group equations in a two- valley system with valley-splitting and intervalley scat- 
tering are derived in the presence of spin-splitting induced by a parallel magnetic field. The relevant 
amplitudes in different regimes set by the relative strengths of the spin and valley splittings and the 
intervalley scattering rate are identified. The range of applicability of the standard formula for the 
magnetoconductance is discussed. 
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I. INTRODUCTION 

In two dimensions, an in-plane magnetic field, -Bj| , 
conveniently couples only to the spin-degrees of free- 
dom leading to spin-splitting of the electron bands. 
Electron- electron (e-e) interactions between the differ- 
ent spin-bands gives rise to a finite magnetoconduc- 
tance, cr(i?|| , T), and hence measurement of cr(S|[ , T) pro- 
vides a simple and accurate way of determining the ef- 
fective spin-related interaction strength^ '^. In a disor- 
dered two-dimensional electron gas (2DEG), the trans- 
port properties at low temperatures, kgT < h/r, are gov- 
erned by singular diffusive particle- hole propagators ^, 
'D{q,uj) = l/{Doq'^ +ui). (Here Dq is the diffusion con- 
stant proportional to the elastic scattering rate l/r.) The 
in-plane magnetic field suppresses the singularity in the 
spin-triplet channels (Sz — ±1) with opposite spin pro- 
jections by introducing a gap proportional to the Zeeman 
energy scale A^ — gfisB^^. This suppression gives rise to 
a negative magnetconductance, Acr(i?j| , T) = cr(i?|| , T) — 
ct(0,T) < in the weak-field limit '*'^ A^ < ksT. In 
the high-field limit, A^ ^ ksT, the spin-bands are well 
split and the transport is governed entirely by the Sz = 
channels which are insensitive to spin-splitting ^^^ . In a 
multivalley system, valley splitting, A„, and intervalley 
scattering A_l = fi/T±, introduces additional gaps ^'^° 
further suppressing cr(_B||,r). The interplay of spin and 
valley degrees of freedom in the presence of finite A„ and 
A_L are studied in this paper. 

It is now well understood that the singular nature of 
I?(g,a;) leads to a strong enhancement of the e-e scat- 
tering amplitudes at low energies ^^. In two dimensions, 
renormalization group (RG) theory applied to a weakly 
disordered system has been extremely successful at cap- 
turing this scale dependence to all orders in the e-e scat- 
tering amplitudes ®'^^. Strong e-e scattering and energy 
renormalization effects, where the latter takes into ac- 
count the renormalization of the Stoner enhancement fac- 
tor, were incorporated into Acr(i3|[ , T) in Refs. 13 and 14. 
They are generalized here to include the effects of A^, and 
A j^ in a two- valley system. [It is assumed that Aj^ < A^, , 
which appears to be the case in high mobility silicon in- 
version layers i^'^^ j ^^j^^ relevant interaction amplitudes 
are identified and the corresponding scaling equations are 
determined in low, A^ < Aj^, and intermediate-fields. 



^± ^ ^z < ^v The high-field hmit A^ > A,„, has 
been studied in great detail in Ref. 17. [Note that for 
ksT < Az, A„, the results in the high and intermediate 
field regimes are the same provided the spin and valley 
indices are interchanged.] 

The RG equations for finite A^ and Aj^ with A^ = 
are detailed in Ref. 18. For finite A^, it is con- 
venient to work in the (spin-singlet/triplets)(8)(valley- 
singlet/triplets) basis so that the spin and valley degrees 
of freedom are treated on equal footing in the particle- 
hole channels. The diffusion propagators and the e-e 
scattering amplitudes in this basis are discussed below. 

Diffusion modes: For A^ — 0, it was sufficient to label 
the modes in terms of the valley states ^'□.(g, w), where 
a = ± and _L. (See Ref. 18 for further details.) a = -I- 
refers to the valley-singlet channel which is gapless, and 
a — — and _L refer to the gapped valley-triplet channels 
with gaps proportional to A_l and A^, -I- A j_ , respectively. 
Since I?_ involves scattering only within the same valley, 
it is insensitive to the splitting A„. It, however, devel- 
ops a gap Aj^ when intervalley scattering mixes the two 
valleys. All modes are effectively gapless at high tem- 
peratures, kgT > A„,A_L. [Temperature units Ty,T± 
and Tz will be used interchangeably in the following to 
represent the scales A^,, A_l and A^, respectively] 

For A^ ^ 0, the spin-triplet channels Sz = ±1 de- 
velop spin-gaps proportional to A^. Hence, it is conve- 
nient to label the propagators as Vta and V^^ , where the 
subscript t corresponds to the spin-triplet channels with 
Sz = ±1, and s labels the Sz = channels, with the 
singlet and the triplet Sz = channels labeled by cr = ±. 

e-e scattering amplitudes: In a single valley system, 
the e-e scattering amplitudes are uniquely described by 
the spin texture of the scattering channel. The ampli- 
tudes Ts and Ft are used to describe the scattering of 
particle-hole pairs in the spin-singlet and triplet chan- 
nels, respectively. They are related to the standard static 
Fermi-liquid amplitudes Fi and r2 as F^. = Fi — F2/2 and 
Ff = — F2/2. These definitions are easily extended to ^^, 
TL = ^2a - 4Tia and Tta = F2Q, where a = ±,_L, 
and (7 — ±. [For notational convenience, the amplitudes 
Fj^ are defined with a factor of -4.] Note that the inter- 
valley scattering amplitudes F^_|_ are generally negligibly 
small in a clean system because the Coulomb scatter- 
ing involving large momentum Qq perpendicular to the 



2D plane is suppressed when the width of the inversion 
layer is many times larger than the lattice spacing, hence 
T'^j_ = rtj_ = T2j_ ■ Together, the total number of ampli- 
tudes equal r-„{4} + TtaH} + r-^{4} + rt±{4} = {16}; 
the number in brackets denotes the number of channels. 
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FIG. 1: Schematics showing the classification of the e-e scat- 
tering amplitudes as a function of temperature T in the pres- 
ence of valley splitting r„ and intervalley scattering Tj_ . The 
relevant amplitudes are marked by solid lines with the degen- 
erate amplitudes grouped together, while dashed lines mark 
the irrelevant amplitudes. The spin and valley structure of 
the amplitudes are described in the main text. 

In the high temperature limit, T > T^,Tz, the am- 
plitudes rj^, except for FJ^, are identically zero. The 
r^_i^ amplitude, which involves scattering in the spin and 
valley singlet channels, (spin-singlet)(g)(valley-singlet), is 
special in that it combines with the long-ranged part of 
the Coulomb interaction to produce a universal ampli- 
tude ^. [Details are given below in section II A.] Hence, 
all 15 of the 16 amplitudes are equal and evolve as r2. 
They are shown grouped together when T > T!u in Fig. 1. 

When Tl ^T <Ty, the 'D± propagators are gapped, 
the corrections to T± are therefore non-singular and 
hence irrelevant. On the other hand, the TJ_ amplitude 
in the (spin-singlet) €5 (valley-triplet) channel, which van- 
ishes at high temperatures, was shown in Ref. 18 to be 
generated under the RG transformations when T < Ty. 
(To emphasize that rj'_ arises as an independent scaling 
variable only when the valley sub-bands are split, it is 
designated as Ty = r^_.) This is a generic feature of 
mutliband systems with sub-band splittings, it was first 
discussed in Ref. 17 in the opposite case Ty < T < T^ 
where the relevant amplitude is T^,^, with the spin and 
valley indices interchanged. 

The splitting of the 15 amplitudes below Ty are shown 
schematically in Fig. 1. The solid horizontal lines mark 
the relevant amplitudes while the dashed lines mark the 
irrelevant ones. The degenerate amplitudes under the RG 
flow are grouped together with the degeneracy indicated 
in curly brackets. At the lowest temperature T < Tl, 
when the two valleys are strongly mixed, only the valley- 
singlet propagator T)^ remains gapless. Hence, only the 
amplitudes in the valley-singlet channel F-|_ survives. 

Clearly, the number of relevant e-e scattering ampli- 
tudes in a multiband system at a given scale is sensitive 
to the splitting and the interband scattering rates within 
the bands. Calculations in different temperature regimes 



in the presence of spin-splitting are detailed below. [Since 
it is seen experimentally that the phase breaking rate 
saturates at low temperatures for low electron densi- 
ties ^^, where the results obtained in this paper are most 
relevant, the contribution from the cooperon (particle- 
particle) channel has been suppressed in the calculations.] 
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FIG. 2: Schematics showing the classification of the relevant 
e-e scattering amplitudes for weak spin-splitting, T^ < Tl, T„, 
as function of temperature. 

Fig. 2 shows schematically the effect of a weak mag- 
netic field Tz < T±,Ty on the amplitudes. The spin 
gap suppresses the singular corrections in the spin-triplet 
channels, hence only TJ_^_ (and TJ'^) develops singular 
diffusion corrections. The amplitude is designated as 
F^ = F^i to emphasize that T ^T^. 
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FIG. 3: Schematics showing the classification of the relevant 
e-e scattering amplitudes for intermediate values of the spin- 
splitting, T± ^Tz <Tv, as function of temperature. 

Fig. 3 shows schematically the relevant amplitudes for 
intermediate values of the spin-splitting T± ^ Tz ^ Ty. 
As in Fig. 2, the spin-triplet channels, Ft, are irrelevant 
below Tz due to the gap in the 2?t propagators. As a 
result, the number of relevant amplitudes (not all degen- 
erate) reduces from seven for T >Tz to three for T <Tz. 



II. SCALING EQUATIONS 

The RG equations in each of the temperature intervals 
shown in Figs. 2 and 3 are derived below. The relevant 
equations when spin-splitting can be ignored, T > Tz, 
have been derived in detail in Ref. 18. The logarithmic 
corrections are presented here in sections HA and II B 
after including the spin degrees of freedom explicitly. 



A. T±,n<T<n 

Since the T>± modes are gapped for T <Ty, their con- 
tributions are non-singular and hence dropped. All other 
modes are effectively gapless when T >Tz,T±. The gap- 



less propagators are set equal to V 



s± 



V, 



t± 



V{q,L0) 



l/{Dq^ + zuj), where D is the renormalized diffusion con- 
stant and z parametrizes the relative scaling of the fre- 
quency with respect to the length scale i^'^o. Both D and 
z acquire diffusion corrections in an interacting system. 
(z = 1 for a non-interacting system ''.) 

The nature of the gapless diffusion modes induce the 
following relations on the amplitudes: Tt+ = Tt- = T2, 
ST+^ = Sr+_ and ST^j^ = 0. Since 6T^j^ = 0, the ampli- 
tudes r^j_ — T2 are degenerate. The diffusion corrections 
in terms of these variables take the form ^* 
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TTiy J (27r)2 
1 f Cpq 



-2T2)V\q,Lo)Dq\ (la) 
r++-2r2)P(q,0), (lb) 



t:v J (27r)2 
cPq 



(r+_ + r++)p(9,o) + 8vi/(r2),(ic) 
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i± 2-KV J (27r)2 



r2p(g,o) + 2*(r2). 



(Id) 



The single integral is defined as / = d^q/(27r)2 and the 
double integral as // = / (Pq/{2i:f J duj/{27r). The den- 
sity of states per spin and valley ly — m/2TT. The con- 
tributions of the "ring" diagrams ^^ equals 4'(r2), where 
(see Fig. 5 in 18) 



*(r2 



r2 [r2p2i 



[rlv^] 



-^ II LuT2[TlV'^] -^Tl[T2V^ 



(2) 



As noted already, the relevance of the T'l_ amplitude in 
the temperature range Jj_ < T < T„ is specific to prob- 
lems with split-bands in a multivalley system. Although 
the corrections 5T^_^_ = <5r^_ for T <Ty, their initial val- 
ues are different. The amplitude r^_ = when T }^ Ty, 
while the singlet amplitude FJ, is special as it combines 
with the static limit of the Coulomb interaction, denoted 
here as FJ^_ (it is conventionally denoted simply as Fq in 
a single valley system with degenerate spin bands ^). The 
r]ti_ amplitude appearing in (la)-(lc) are to be replaced 
by its long-ranged value 



F+ 



^1 — -L 04 



F+ 



(3) 



When combined with the F2 amplitude, the long-ranged 
singlet amphtude is given as: Ff ^ = F2 - 4Ff'^. (To 



be consistent with the notations in this paper, T^"- is 
defined with an extra factor of -4.) It is easily verified 
by combining Eqs. (Ib)-(ld) that the singlet combination 
6{z + T^^) = is satisfied at all length scales, provided 
the corrections to the static amplitude STq_^_ — 0. This is 
a well established result, with great importance for the 
general structure of the theory H'l^ 

Having obtained the leading logarithmic corrections, 
the scaling equations are derived to first order in the 
dimensionless resistance p = l/4(27r2j/£)) and to all or- 
ders in the e-e interaction amplitudes by performing the 
ladder summations described in Fig. 6 in Ref. 18. It 
amounts to replacing the static amplitudes F^ by the dy- 
namical amplitudes Ui{q,Lu): 



Ui{q,uj) 



2?j(g,a;) 



''V{q,Lo)' 
where, the propagators Vi are defined as 

1 



Dq^ + {z + Ti)uo 



(4) 



(5) 



The amplitudes F^ represents F2,F^_ and Ff'-'^. Note 
that since the leading logarithmic corrections involve only 
one momentum integration, it generates only one factor 
of 1/D. The corrections are therefore limited to the first 
order in resistance p (disorder). The limitation on the 
number of momentum integrations constraints the num- 
ber of e-e vertices in the skeleton diagrams. The lad- 
der sums extend the skeleton diagrams to all orders in 
Ti without changing the number of momentum integra- 
tions. Note, however, that only those interaction vertices 
involving frequency integrations can be extended to in- 
clude dynamical effects. These amplitudes are enclosed 
in square brackets in Eq. (2). Substituting the F2 ampli- 
tudes in the square brackets with U2 and performing the 
q,aj integrals leads to the very simple expression ^'^^: 



*(F2) 



?)^f- 



1 
TV 



(6) 



The remaining single integrals J cP'q'D{q, 0) involving 
only momentum integrations are easily evaluated to give 



<Pq 

J2n) 



^V{q,0)^2p\og 



1 
TV 



(7) 



The integrals in SD containing lu integrations remain 
to be evaluated. Before the integrals can be done, the F^^_ 
amplitude is replaced with Ff^ following Eq. (3), after 
which the amplitudes Ff ^, F^_ and F2 are rearranged to 
form T^^ and F+_ and F2 and then extended to C/f^, U^_ 
and C/2, respectively. 

When the equation for p is expressed in terms of the 
scaling variables, 72 — F2/2; and 7„ = F^_/z, the equa- 
tions for yO, 72 and 7^ form a closed set of equations in- 
dependent of z. The final RG equations in the range 



T±,Tz ^T <Ty are given below with the scale ^ defined 
to logarithmic accuracy as ^ = log(l/rr). 



dp 

drf2 
d^ 

rflnz 



= p\l-^^^)-mi2)), (8a) 

= f [(l + 72)' + (l + 72)(72-7.)], (8b) 



= §(l+7«)(l-7. 



672), 



2 ^ 



Iv - 672) 



(8c) 
(8d) 



The function $(7) is defined as 



$(7)= (l + -)log(l + 7)-l 



(9) 



As described in Fig. 3, the 15 degenerate amplitudes 
for T > Ty split into six r2 and one r„ amplitude when 
T ^ Ty. [This splitting of the amplitudes is generic to 
multiband systems with subband splittings. The same 
equations are obtained when instead of the valley bands, 
the spin bands are split ^^, i.e., T^ <T < T^] Note that 
7„ coincides with 72 when T ^ Ty. 



B. T,<T<T±,n 

The relevant amplitudes in the presence of strong val- 
ley mixing (T < T±) correspond to scattering in the 
valley-singlet channels, FJ^ and Tt+ ■ Since ST'^_^ ~ 
vanishes in the absence of spin-splitting {T >Tz), it fol- 
lows that the amplitudes F~^ — Ft+ — F2 are all equal 
and satisfy the equations ^^ 



5D 

Id 

5z 
(5F2 



{T+^-T2)V''{q,uj)Dq\ (10a) 

.W i^An,-T.)Viq,,), (10b) 

^/(0r++2?('?,O)+4vI/(F2), (10c) 
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d^q 
Attv J (27r)2 



F22?(g,0) + *(F2). (lOd) 



The coefficient of F2 and the ring diagrams ^(F2) in 
(lOa)-(lOc) are suppressed by a factor two when com- 
pared with (la)-(lc) since they no longer contain a val- 
ley sum. The corrections to FJ'^ in (lOd) already do 
not carry a valley sum, only half the amplitude involving 
the same valley, however, acquires corrections when the 
valley-bands are mixed, which accounts for the overall 
factor of half when compared with (Id). Note that the 
condition 5{z + Ff'^) = is satisfied. Following the pro- 
cedure described in section II A, the RG equations read 



dp 
d^ ^ 


= p2 (1 _ 3<i>(72)) 


(11a) 


d72 
d^ 


= ^(1+72)^ 


(lib) 


d\nz 

dS 


= -f(l-372) 


(lie) 



The function $(7) is defined in Eq. (9). As described in 
Fig. 2, only three of the 15 degenerate amplitudes survive 
when T < T± when spin-splitting can be neglected T > 
Tz- The high field cases are discussed below, i.e., T <Tz. 



C. T±<T<n<T^ 

It should be noted that the results for T <Tz "^T^ is 
equivalent to the situation if the gap scales were reversed, 
i.e., T < Ty < Tz, provided of course the spin and valley 
indices are interchangeable, which is the case when T > 
T±. The RG equations for T^ < T^ are derived in Ref. 17. 
The opposite situation T^ < T^ is derived here. 

When T < Tz, the 'Dt± propagators are gapped, and 
hence the corrections in the 5'^ = ±1 channel are non- 
singular. The corresponding amplitudes Tt± are there- 
fore irrelevant at these temperatures, which reduces the 
number of relevant interaction amplitudes by four. Fur- 
thermore, the amplitude Fj~_^ acquires diffusion correc- 
tions ^^ when T < T^ in the same way that r^_ does 
when T <Ty. Since T <Tz and Ty, the amplitude F^_ 
also acquires logarithmic corrections. As a result, both 
F^j_ are different from F2 when T < T^ , T^ . After includ- 
ing the contributions from FJ^,, the diffusion corrections 
for T < T, take the form: 



SD 

Sz = 

sr2 = 

"^ la — 



ljj[T. rL-r2J p3(g,^)Dg^(12a) 

— ((T. rL-r2) v{q,o), (i2b) 

^ / E rL2?(9,0)+4*(F2), (12c) 

— /F2l?(g,0) + *(F2). (12d) 



The coefficient of F2 and the ring diagrams ^(F2) in 
(12a)-(12c) are suppressed by a factor two when com- 
pared with (la)-(lc) since they no longer contain a spin 
sum. The corrections to FJ^ in (12d) already do not 
carry a spin sum. Only half the amplitude involving the 
same spin, however, acquires corrections when the spin- 
bands are split, which accounts for the overall factor of 
half when compared with (Id). Since <5Fj~^_ = i5Fj~_, the 
amplitudes, after combining with F2, can be grouped to- 
gether as Tj^ = Tz . Extending the singlet amplitude F^_^ 



to include the static long ranged part of the Coulomb 
interactions Ff ^ as discussed in Eq. (3) and using the 
identity r2 - E„,, rf„ - Ea,. ^1J^. the amplitude T^ 
can be eliminated from Eqs. (12a)-(12d) in favor of the 
amplitudes T^'^.Tz and r„ as 



6D 

Id 

5z 



\ JJi^s" + r. + 2T,)V^iq,Lu)Dq^, (13a) 

(rf^^ + r,„ + 2r,)p(q,o), (i3b) 

(13c) 



1 

4:TTl> 

sr„ = (5rf'« = -s 



Combining (12c) and (12d) to give Eq. (13c) is possi- 
ble only because the ^(r2) contribution cancels exactly 
when the sum over opposite spin-projections are sup- 
pressed due to spin-splitting '^ . Also note in (13c), that 
the singlet combination S{z + T^^) = holds explicitly, 
as needed for the consistency of the RG theory^' ^^. 

The RG equations are obtained by evaluating the inte- 
grals after extending the static amplitudes by the dynam- 
ical amplitudes Ui defined in Eq. (4). The RG equations 
for T^<T<T^<Ty are 



P'il 



$(7.) - 2<i>(7,)) , 
(l + 7,)(l-7,„-27,) 
(l+7„)(l-7. -27.), 



-f (1-7. -27.) 



(14a) 
(14b) 
(14c) 
(14d) 



dp 
d^ 

djz 

di 
di 

dim 

As described in Fig. 3, the four Tt± 
suppressed when T < T^, leaving two 
which evolve away from F2. Note that 7^ 
Tz , while jy w 72 when T ^ Ty. [The RG equations when 
spin-splitting is large, T < Ty ^ T^, take the same form 
as Eqs. (14a)-(14d) provided the spin and valley indices 
are interchanged; see Ref. 17 for details.] 

D. T<T^,T,,n 

The two valleys are strongly mixed when T < T±, 
leaving only the valley-singlet propagators V^^ gapless. 
Hence, only F7+ = F^ and F^^, survive at low tempera- 
tures. The corresponding diffusion corrections read 



amplitudes are 
Tz amplitudes, 
w 72 when T w 



W 

Sz 
5F2 



sn+ 



III (rr+ + r^+ - ^r2) v\q,uj)Dq\ (isa) 
-l^ J (rr+ + r^+ - ^rs) v{q, o), (isb) 
;^/(rr+ + r++)i?(<z,o) + 2*(r2), (i5c) 

^jT2V{q,0) + ^^ir2) (15d) 



All terms involving F2 amplitudes are suppressed by 
a factor of two in (15a)-(15c) compared to (12a)-(12c) 
due to the suppression of the Ft_ amplitudes, which are 
irrelevant when T < T± . The equations can be simplified 
in terms of the amplitudes T^^ and F^ as 



6D 


- 1 /■/ 


D 


- ^77 


5z -- 


1 


5Vz -- 


= 5T':" 



{T'^^'' + Tz)VHq,u:)Dq\ 


(16a) 


f{r^^ + Tz)V{q,0), 


(16b) 


= -6z. 


(16c) 



Note again that the condition 5{z + F^^) = is sat- 
isfied. Following the procedure followed in the previous 
sections, the RG equations for T <T±,Tz,Ty are 



dp 
d^ 


= p^l~^Jz)), 


(17a) 


diz 
d^ 


= f(l + 7.)(l-7.), 


(17b) 


dlnz 
d£ 


^-f(i-..)^ 


(17c) 



These equations coincide with the results obtained in 
the case of a single valley with spin-splitting studied in 
Ref. 20. Strong intervalley scattering for T < T± mixes 
the two valleys to effectively produce a single valley. 



III. CONCLUSIONS 

The derivation of the scaling equations in Sec. II were 
carried out keeping only the gapless valley and spin chan- 
nels in each temperature interval. The scale dependence 
of the dimensionless resistance p = (e^/7r/i)i?n, where 
R\j is the sheet resistance, is then obtained by integrat- 
ing the self consistent set of scaling equations separately 
in each temperature interval and matching the values of 
the amplitudes and resistance at the boundaries of each 
interval. Since the intervals are sensitive to the value of 
Tz, one obtains in this way p{B\\,T) whose inverse gives 
cr(i?||,T) = l/p(i3||,T). This method is, however, not ac- 
curate as the crossover regions have finite contributions 
from the gapped channels and are hence non-universal. 

The case of weak spin-splitting T^ <T can be treated 
fairly accurately, however. In this case the sensitivity to 
Bu arises only from the presence of a weak spin-gap in 
the triplet channels below the scale set by T. Hence, 
subtracting a{0, T) from <y{B\\ , T) captures only the con- 
tributions originating from the suppression of the triplet 
channels. The explicit form Aa{B\\,T) for the single val- 
ley case was derived in Refs. 5,14 in the limit T^ ^ T. 
It is straightforward to extend the results to include the 
valley degrees of freedom by keeping track of the number 
of F( amplitudes in a given temperature interval as ^ 



Aa(i3,i,r) 



-0.091--iVt72(72 + l)(r,/T)2. (18) 
27r/i 



The number of relevant Fj amplitudes, Nt can be read off 
directly from Figs. 1, 2 and 3 as TV* = 2 for T;, < T < Ti, 
TVi = 4 for T; < T < T^ and Nt ^ 8 ior T > T^,T^. 

To summarize, RG equations in the presence of spin- 
splitting induced by a parallel magnetic field have been 
obtained in a two valley system in the valley-split and 
strong intervalley scattering regimes. The behavior of 
the magnetoconductance in the limit of weak magnetic 
field are discussed. 
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